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CUBATURE  FORMULAS  AND  RELATED  QUESTIONS1 


V.  N.  Temlyakov 

Department  of  Mathematics,  University  of  South  Carolina,  Columbia,  SC  29208 

E-mail:  temlyak@math.sc.edu 


Abstract.  The  main  goal  of  this  paper  is  to  demonstrate  connections  between  the 
following  three  big  areas  of  research:  the  theory  of  cubature  formulas  (numerical 
intergation),  the  discrepancy  theory,  and  nonlinear  approximation.  In  Section  1  we 
discuss  a  relation  between  results  on  cubature  formulas  and  on  discrepancy.  In  par¬ 
ticular,  we  show  how  standard  in  the  theory  of  cubature  formulas  settings  can  be 
translated  into  the  discrepancy  problem  and  into  a  natural  generalization  of  the  dis¬ 
crepancy  problem.  This  leads  to  a  concept  of  the  r-discrepancy.  In  Section  2  we 
present  results  on  a  relation  between  construction  of  an  optimal  cubature  formula 
with  m  knots  for  a  given  function  class  and  best  nonlinear  m- term  approximation 
of  a  special  function  determined  by  the  function  class.  The  nonlinear  m- term  ap¬ 
proximation  is  taken  with  regard  to  a  redundant  dictionary  also  determined  by  the 
function  class.  Sections  3  and  4  contain  some  known  results  on  the  lower  and  the 
upper  estimates  of  errors  of  optimal  cubature  formulas  for  the  class  of  functions  with 
bounded  mixed  derivative.  One  of  the  important  messages  of  this  paper  is  that  the 
theory  of  discrepancy  is  closely  connected  with  the  theory  of  cubature  formulas  for 
the  classes  of  functions  with  bounded  mixed  derivative.  We  have  included  in  the 
paper  both  new  results  and  known  results.  We  included  some  of  known  results  with 
their  proofs  for  the  following  two  reasons.  First  of  all  we  want  to  make  the  pa¬ 
per  selfcontained  (within  reasonable  limits).  Secondly,  we  selected  the  proofs  which 
demonstrate  different  methods  and  are  not  very  much  technically  involved.  Section 
5  contains  historical  notes  on  discrepancy  and  cubature  formulas,  some  further  com¬ 
ments  and  remarks.  Historical  remarks  on  nonlinear  approximation  are  included  in 
Section  2.  We  want  to  point  out  that  this  paper  is  not  a  survey  in  any  of  the  above 
mentioned  areas.  We  did  not  even  try  to  provide  a  complete  list  of  results  in  those 
areas.  We  rather  wanted  to  highlight  the  most  typical  results  in  cubature  formulas 
(Sections  3  and  4)  and  show  their  relation  to  the  discrepancy  theory. 


1.  Cubature  formulas  and  discrepancy 

Numerical  integration  seeks  good  ways  of  approximating  an  integral 


f(x)d[i 


1This  research  was  supported  by  the  National  Science  Foundation  Grant  DMS  9970326  and 
by  ONR  Grant  N00014-96-1-1003 
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by  an  expression  of  the  form 

m 

Am(/,0  :=!>;/(?),  t  =  e'ea  j  =  (1.1) 

3  = 1 

It  is  clear  that  we  must  assume  that  /  is  integrable  and  defined  at  the  points 
The  expression  (1.1)  is  called  a  cubature  formula  (A,  £)  (if  Q  C  Kd, 
d  >  2)  or  a  quadrature  formula  (A,  £)  (if  Oct)  with  knots  £  =  (£x, . . .  ,  £m)  and 
weights  A  =  (Ai, . . . ,  Am).  For  a  function  class  IT  we  introduce  a  concept  of  error 
of  the  cubature  formula  Am(-,£)  by 

Am(IF,0  :=  sup  |  [  fdp  -  Am(/,0|-  (1-2) 

/gw  Jn 

In  order  to  orient  the  reader  we  will  begin  with  the  case  of  univariate  periodic 
functions.  Let  for  r  >  0 

OO 

Fr(x)  :=  1  +  2  E]  k  r  cos (kx  —  m/2)  (1.3) 

k= 1 


and 

W;  :={/:/  =  M\P<1}  (1.4) 

where  *  means  convolution  and  ||  •  ||p  is  the  standard  Lp-norm.  It  is  well  known  that 
for  r  >  1/p  the  class  W£  is  embedded  into  the  space  of  continuous  functions  C'(T). 
In  a  particular  case  of  W/  we  also  have  embedding  into  C'(T).  From  the  definitions 
(1.1),  (1.2),  and  (1.4)  we  see  that  for  the  normalized  measure  dp  =  ^ dx 


a™(it;,o 


sup  [  Fr{x -y)dp-S"\jFr(^  -y))p{y)dy\ 

M„<i  l7T  J t  Jt 


=  111  -  E  XiFr(?  ~  Oil ✓  >  p’  ■■=  (1-5) 

Thus  the  quality  of  the  quadrature  formula  Am(-,£)  for  the  function  class  W£  is 
controlled  by  the  quality  of  Am(-,£)  for  the  representing  kernel  Fr(x  —  y ).  In  the 
particular  case  of  W\  we  have 

m 

Am(IT11,0  =  max  1 1  -  Va^Ti^  -  y)  |.  (1.6) 

y  ft 


In  this  case  the  function 

OO 

F\  (s)  =  1  +  2  sin  kx  =  1  +  S(x) 

k=  1 


has  a  simple  form:  S(x)  =  0  for  x  =  In  and  S(x)  =  tt  —  x  for  x  G  (0,  2tt). 
This  allows  us  to  associate  the  quantity  Arn  ( W/ ,  £ )  with  the  one  that  has  simple 
geometrical  interpretation.  Denote  by  %  the  class  of  all  characteristic  functions 
X[o ,a]ix),  a  £  [0,  2tt).  Then  we  have  the  following  property. 

2 


Proposition  1.1.  There  exist  two  positive  absolute  constants  C\  and  C2  such  that 
for  any  Am(-,£)  with  a  property  'fT  .  A j  =  1  we  have 

C,iAm(X,0  <  A m(Wl,0  <  C2 Am(x,0-  (1-7) 


Proof.  We  have  for  any  a  G  [0,  27t) 


X[o,„]M  =  ^-(<J  +  -FiM  -  F1(x-a)). 


(1.8) 


Thus  using  (1.6)  we  get 


Am(x,0<-Am(^i1,0, 


7T 


which  proves  the  left  inequality  in  (1.7). 

Let  us  prove  the  right  inequality  in  (1.7).  Denote  e  :=  Am(x, £).  Then  by  (1.8) 
we  get  for  any  a  G  [0,  2n) 


—27 re  <  (Fi(x)  —  Fi(x  —  a))dn  —  Am(Fi(x)  —  Fi(x  —  a),  £)  <  2ire.  (1.9) 

Jt 

Intergating  these  inequalities  against  a  over  T  we  get 

|  [  F1(x)d/j,  -  Am(F\,£)|  <  2?re. 

J  T 

This  inequality  combined  with  (1.9)  implies 

Am{WlO  <47rAm(X,0- 


We  proceed  to  the  multivariate  case.  For  a;  =  (aq, . . . ,  Xd)  denote 

d 

Fr{x)  :=  Fr(xj) 

3= 1 

and 

MWp  :={/:/  =  <^  *  Fr,  \\<p\\p  <  1}. 

For  /  G  MWp  we  will  denote  /(r)  :=  99  where  99  is  such  that  /  =  9 p  *  Fr.  The 
letter  M  in  the  notation  MWf  stands  for  ’’mixed”,  because  in  the  case  of  integer  r 
the  class  MWf  is  very  close  to  the  class  of  functions  /,  satisfying  ||/(r,-"’r)  ||p  <  1, 
where  f(r’~-’r'>  is  the  mixed  derivative  of  /  of  order  rd.  A  multivariate  analog  of 
the  class  x  is  the  class 

d 

Xd  ■=  {X[o,a](aO  :=  IJ  X[o,aj](xj),  aj  G  [0,  2tt),  j  =  l,...,d}. 

7  =  1 

Similarly  to  the  univariate  case  one  obtains  analogs  of  (1.5)  and  (1.6) 

m 

A =  111  -  J]  A jFr(?  -  OIU;  (1.10) 

3  = 1 

m 

A  m{MWlO=  max  1 1  -  V  A  j  F1  (£>'  -  y)  |  • 

W  i=i 

We  prove  a  multivariate  analog  of  Proposition  1.1. 
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(1.11) 


Proposition  1.2.  There  exist  two  positive  constants  C\(d)  and  C2(d)  such  that 
for  any  Am(-,£)  with  a  property  fT  .  A j  =  1  we  have 

CMAmix^O  <  A m{MWle>  <  C2(d) Am(xdiO-  (1-12) 


Proof.  First  we  prove  the  left  inequality.  Denote  e  :=  Am(MWf,  £).  Then  by  (1.11) 
we  have 

m 

-e<l-J2^Fi(^~y)<e,  yeTd. 

n=i 

Take  any  subset  e  C  [1,  d]  and  integrate  the  above  inequality  against  y3 .  j  e,  over 
Td-Iel .  We  get 

m 

Thus  for  any  e  C  [1,  d],  e  ^  0,  and  any  yj,  j  G  e,  we  have 

m 

i  e  aaii  F^i  -  -  n  -  n))  i  s  2f. 

M=1  j Ge  j£e 

This  inequality  and  the  representation  (1.8)  imply  the  left  inequality  in  (1.12). 

We  prove  the  right  inequality  in  (1.12)  by  induction.  In  the  case  d  =  1  it 
follows  from  Proposition  1.1.  Assume  we  have  it  for  all  dimensions  d!  <  d.  Denote 
e  :=  Am(xd,£).  For  a  subset  e  C  [1  ,d]  we  denote 

«e):=K1 * 3(e),...,f"(e)),  E'lelet1'1,  ("(e) j  =  j  e  e. 

Consider  the  class  x^(e)  that  is  the  class  ylel  corresponding  to  the  coordinated  Xj 
with  j  G  e.  It  is  clear  that  we  have 

A m(x|e|(e),^(e))  <  e. 

Therefore,  by  the  induction  assumption  and  by  (1.11)  we  have  for  all  e  c  [l,d], 
\e\  <  d ,  that 


I1 -E^IPi «?-»)!  SC(d)e,  s,  eT,  jet  (1.13) 

M=1  j£e 

It  is  easy  to  see  that  for  functions  of  the  form 

cl 

fa,y(x )  :=  II  (X[0>t»3-]  (a=j  -  %•)  -  7^) 

3  = 1 
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where  X[o ,aj](xj  ~  Uj)  is  27t  periodic  we  have 


\Am(fa,y,0\<C(d)e. 

Using  this  estimate  and  the  representation  (1.8)  we  obtain 

d 

lAm(JJ(^i(a:i  ~Vj)  ~F1(xj  -  yj  -%))!  <  C(d)e 
j  i 

what  means 

m  d 

-C(d)e  <  K  n  F^J  -  Vi) 

n=l  j= 1 
m 

+  E  EAdIFiK?-w)n<-FiK?-»-%))SCMe. 

e:|e|<d/x=l  j£e  j(^e 

Integrating  against  a  over  Td  we  get  from  here 

m  d  m 

i  e  a-  n  -  vi)  +  e  (-1)"-1*1  e  n  -  n)\  <  c(d)e. 

M=1  i=1  e:\e\ <d  /i=l  3^e 

Using  (1.13)  and  the  identity 

E  (-i)d-|-i  =  -i 

e:\e\ <d 


we  obtain 

m  d 

E^n^Kl-W)-1!  <C(d)e. 

11=1  j  =  1 

We  complete  the  proof  by  applying  (1.11). 

The  classical  definition  of  discrepancy  (in  the  convenient  for  us  form)  of  a  set  X 
of  points  x1, . . . ,  xm  C  [0,  l]d  is  as  follows 


D(X,m1d)00 


max 

a<E[0,l]d 


d  1  m 

riai  -  Mi- 


it  is  clear  that 


D(X,m,  d)oo  =  Am(xd,27iW)  with  Ai  =  •  •  •  =  Xm  =  1/m. 

Thus  by  Proposition  1.2  the  classical  concept  of  discrepancy  is  directly  related  to 
the  efficiency  of  the  corresponding  cubature  formula  for  a  special  function  class 
MWl .  It  is  well  known  that  the  Wj  is  very  close  to  the  class  of  functions  of 
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bounded  variation  and  the  MWj  is  very  close  to  the  class  of  functions  with  bounded 
variation  in  the  sense  of  Hardy- Vitali.  In  the  beginning  of  20th  century  D.  Vitali  and 
G.  Hardy  generalized  the  definition  of  variation  to  the  multivariate  case.  Roughly 
speaking,  in  the  one-dimensional  case  the  condition  that  /  be  of  bounded  variation 
is  close  to  the  condition  1 1 /' 1 1 1  <  oo.  In  the  multidimensional  case  the  condition 
that  a  function  have  bounded  variation  in  the  sense  of  Hardy-Vitali  is  close  to  that 
requiring  ||i  <  oo,  where  /C1’---’1)  is  a  mixed  derivative  of  /. 

In  the  thirties  in  connection  with  applications  in  mathematical  physics,  S.L. 
Sobolev  introduced  the  classes  of  functions  by  imposing  the  following  restrictions 

||j(nlv..,nd)||p  <  x  (1.14) 

for  all  n  =  (ni, . . . ,  n^)  such  that  n\  +  •  •  •  +  rid  <  R.  These  classes  appeared  as 
natural  ways  to  measure  smoothness  in  many  multivariate  problems  including  the 
numerical  integration.  It  was  established  that  for  Sobolev  classes  the  optimal  error 
of  numerical  integration  by  formulas  with  m  knots  is  of  order  m~Rld .  Assume 
now  for  the  sake  of  simplicity  (to  avoid  fractional  defferentiation)  that  R  =  rd , 
r  natural  number.  At  the  end  of  fifties,  N.M.  Korobov  discovered  the  following 
phenomenon.  Let  us  consider  the  class  of  functions  which  satisfy  (1.14)  for  all  n 
such  that  rij  <  r,  j  =  1, . . . ,  d  (compare  to  the  above  classes  MW J).  It  is  clear 
that  this  new  class  (class  of  functions  with  dominating  mixed  derivative)  is  wider 
then  the  Sobolev  class  with  R  =  rd.  For  example,  all  functions  of  the  form 

d 

f(x)  =  '[[fj{xj),  ||/jr)||p<l, 

3  =  1 


are  in  this  class,  while  not  necessarily  in  the  Sobolev  class  (it  would  require,  roughly, 
ll/j^llp  <  1)-  Korobov  constructed  a  cubature  formula  with  m  knots  which  guar¬ 
anteed  the  accuracy  of  numerical  integration  for  this  class  of  order  m_r(logm)rd, 
i.e.  almost  the  same  accuracy  that  we  had  for  the  Sobolev  class.  Korobov’s  discov¬ 
ery  pointed  out  the  importance  of  the  classes  of  functions  with  dominating  mixed 
derivative  in  fields  such  as  approximation  theory  and  numerical  analysis.  The  con¬ 
venient  for  us  definition  of  these  classes  (classes  of  functions  with  bounded  mixed 
derivative)  is  given  above  (see  the  definition  of  MW J). 

In  addition  to  the  classes  of  27r-periodic  functions  it  will  be  convenient  for  us  to 
consider  the  classes  of  nonperiodic  functions  defined  on  Qd  =  [0,  l]d. 

Let  r  be  a  natural  number  and  MWp(0,d ),  1  <  p  <  oo  denote  the  closure  in 
the  uniform  metric  of  the  set  of  rd-times  continuously  differentiable  functions  f(x ) 
such  that 


MW;(Qd) 


E 


Qn1  +  - 


•+nd 


f 


dx ^  . . .  dxndd 


<  1. 


j  =  l,...,d 


llsllp 
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where 


It  will  be  convenient  for  us  to  consider  the  subclass  MWf(Qd)  of  the  class 
MWp(Qd )  consisting  of  the  functions  f(x)  representable  in  the  form 


where 


f(x)=[  Br(t,x)(p(t)dt ,  \\(p\\p  <  1, 

J  fid 


d 

Br(t,  x)  =  JJ((r-  1)!)_1(^  -Xj)7^1 
i= 1 

t,  x  G  0^,  (a)+  =  max(a,  0). 


In  the  connection  with  the  definition  of  the  class  MWp{Q,d)  we  remark  here  that 
for  the  error  of  the  cubature  formula  (A,  £)  with  weights  A  =  (Ai,...,Am)  and 
knots  £  =  (£1, . . . ,  £m)  the  following  relation  holds.  Let 


A m(/,£)  -  [  f(x)dx 
J  fid 


then  similarly  to  (1.5)  and  (1.10)  one  obtains  (p1  =  p/(p  —  1)) 


Am{Mw;(nd),c)  =  sup  Rm(A,£,f)  = 

f<EMW£(fld) 

m  d 

M=1  3= 1 

The  quantity  Hr(£,  A,  m,  d)q  in  the  case  r  =  1,  A  =  (1/m, . . . ,  1/m)  is  the  clas¬ 
sical  discrepancy  of  the  set  of  points  {£M}.  In  the  case  A  =  (1/m, . . . ,  1/m)  we 
denote  Dr(£,m,d)q  :=  Dr(£,  (1/m, . . . ,  1/m),  m,  d)q  and  call  it  the  r-discrepancy. 
Thus  the  quantity  _Dr(£,  A,  m,  d)q  defined  in  (1.15)  is  a  natural  generalization  of 
the  concept  of  discrepancy.  This  generalization  contains  two  ingredients:  general 
weights  A  instead  of  the  special  case  of  equal  weights  (1/m, . . . ,  1/m)  and  any  nat¬ 
ural  number  r  instead  of  r  =  1.  We  note  that  in  approximation  theory  we  usually 
study  the  whole  scale  of  smoothness  classes  rather  than  an  individual  smoothness 
class.  The  above  generalization  of  discrepancy  for  arbitrary  positive  integer  r  al¬ 
lows  us  to  study  a  question  of  how  does  smoothness  r  affect  the  rate  of  decay  of 
generalized  discrepancy. 

Remark  1.1.  In  the  case  of  natural  r  the  class  MWf  turns  into  the  subclass  of  the 
class  MWp(Qd)B  :=  {/  :  f  / B  G  MW£ (fid)},  after  the  linear  change  of  variables 

Xj  =  —7 r  +  2irtj,  j  =  1, . . . ,  d. 


=  A-(£,A,m,dy.  (1.15) 


We  are  interested  in  dependence  on  m  of  the  quantities 


5m(W)=  inf  A  m(W,0 


for  the  classes  W  defined  above. 

It  will  be  convenient  for  us  to  use  the  following  notations.  For  two  nonnegative 
sequences  {am}^=1  and  1  we  write  am  bm  or  brn  am  if  there  exists 

a  positive  C  independent  of  m  such  that  am  <  Cbm ,  m  =  1,2, ... .  In  the  case 
am  -C  bm  and  am  bm  we  write  am  x  bm. 

Remark  1.1  shows  that 

Sm(MW; )  «  5m(MW;(Qd)).  (1.16) 


0 

Let  MWp(0,d)  denote  the  subset  of  functions  in  MWf(Lld )  which  is  the  closure 
in  the  uniform  metric  of  the  set  of  functions  /  which  have  the  following  property: 
f(x)  is  rd  times  continuously  differentiable  and  suppf(x)  C  f Id- 

Theorem  1.1.  Let  1  <  p  <  oo.  Then 

5m(MWrp(nd ))  X  5m(MW;{Qd))  X  5m(MW;(Qd)).  (1.17) 


Proof.  Let  A  and  f  be  given.  We  will  prove  that 


sup 

g£MWp(fld) 


»  sm(Mw;(sid)). 


(1.18) 


Suppose  an  infinitely  differentiable  function  tf(x)  is  such  that  *f(x)  =  0  for  x  <  0, 

■0( x )  =  1  for  x  >  1  and  ip(x)  strictly  increases  on  [0, 1].  For  the  cubature  formula 

o 

(A,  £),  defined  on  the  class  MWp(Lld),  we  define  the  cubature  formula  (A ' ,rj),  whose 
error  will  be  investigated  for  the  class  MWp(Lld )  as  folows: 

j  =  l,...,d; 

d 

^  =  ^11^)’  /i  =  1,  .  .  .  ,  TO. 

3  = 1 

Then  for  the  functions  /  and  g  related  as 


9(t) 


d 


/(V’(ti),...  ,^(td))  Yl'if'itj) 

3  = 1 


(1.19) 


we  have 


f(x)  dx=  g(t )  dt, 


Qd 


Qd 


E  w)  =  E  vw*)- 

!d=l  g,=  l 
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It  remains  to  check  that  there  exist  a  number  S  >  0  which  does  not  depend  on  m 

such  that  5g  G  MWp(Qd)  provided  /  G  MW f  (fid)- 

In  fact,  differentiating  (1.19)  we  see  that  the  expression  for  g ^  (t),  s  =  (si, . . . ,  Sd), 
0  <  Sj  <  r,  j  =  1, . . . ,  d  will  contain  terms  of  the  form 

d  i 

u{t,k) = /(fc) (vm) nn^wr.  ^ = 

3=  1  * 

k  (^i)  ■  ■  ■  j  k ci ) ,  0 . .  kj  . Sj ,  ^  '  I'jTfij  Sj  T  1,  j  1)  ■  ■  ■ )  d 

i  (1.20) 


the  number  of  which  depends  on  the  vector  s. 
It  is  obvious  that  in  the  case  p  =  oo  we  have 


u{t,k) ||oo  <  C(r)\\f\\Mw^{nd)- 


To  estimate  ||oj(-,  k)  ||p,  1  <  p  <  oo  we  use  the  following  simple  lemma. 

Lemma  1.1.  Suppose  that  f  G  MWf(Qd)  is  rd  times  continuously  differentiable, 
and  the  vector  k  G  is  such  that  kj  =  r  for  j  G  e\  and  kj  <  r  —  1  for  j  E  e 2  = 
[1,  d]  \ei.  Then 


sup 

Xj,j€e2 


J  | /<*"• 

[0,l]leil 


Wdl  dx^  -  C'CP*  rf)  1 1  / 1 wj  (n„ 

i£e  1 


1  <  p  <  00. 


Proof.  We  will  first  prove  the  following  statement.  Let  /  be  such  that  /,  are 
continuous.  Then 


sup  f  \f(x)\p{Y[dxi)  <2p(\\f\\P  + 

Xj  J  . 

[0)i]d-i 


(1.21) 


Indeed,  there  is  an  a  G  [0, 1]  such  that 


/ 


[0 , 1] d  1 


|/(Z1,...  ,Xj 


i,a,xj+1,...  ,xd)\p{Y[dxi)  <  \\f\\p. 
i¥=j 


We  represent  the  function  f(x)  in  the  form 


f  ( x )  f  (x  1 ,  •  •  •  ,  Xj — 1 ,  a,  Xj  —  1 , . . .  ,  ^T(i) T 


1 ,  ’p/  —  1  j  ■  ■  •  5  ^d)  du. 
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Then,  for  any  Xj  G  [0, 1],  we  have 


This  proves  (1.21). 

Applying  successively  relation  (1.21),  we  obtain  the  statement  of  the  lemma. 

We  return  to  estimating  ||cl»(-,  A:)||p.  By  Lemma  1.1  and  by  the  uniform  bound¬ 
edness  of  the  functions  |'0^(t7-)|  <  C(r),  l  <  r  +  1,  we  get  for  k  such  that  k3  <  r, 
j  =  1, . . . ,  d  that 

\\u{-,k)\\p  <  |M-,&)||oo  «  ||/(fcl’"',fed)||oo  <  \\f\\MW£{nd)- 

Thus  it  remains  to  consider  only  those  k  for  which  there  is  a  j  such  that  kj  =  r. 
Then,  with  respect  to  this  variable,  {'^'(tj))r+1  participates  as  an  additional  cofactor 
in  expression  (1.20).  Taking  into  account  that 

M’'te))’,<r+1)  <C(p,rW(tj), 

we  obtain 

l 

J \co(t,  k)\p  dtj  < 
o 
i 

<C{p,r)( /|/(^i),..  ,il>(tj-1),xj,il>(tj+ 1),...  ,il>(td))\p  dxj)x 
o 

xnnw(0(«)mtP- 

v^j  i 

Reasoning  in  this  way  for  all  j  such  that  kj  =  r  and  applying  Lemma  1.1,  we  find 
that 

M-,fc)||p  <  C{p,r,d)\\f\\Mw-(Qd) 

o 

for  all  k.  This  implies  that  c(p,r,d)g  G  MWrp{Q,d)  with  some  positive  c(p,r,d). 
The  arguments  we  presented  show  that 

P  m 

sup  I  j  g(t )  dt-^2  A pg(^)  \  > 

g£MWp(fld)  Qd  M=1 


>  c(p,  r,  d)  sup 

f€MW;(nd) 


/  f(x)  dx 
Qd 


Eb/w*)  i- 

n=l 


(1.22) 


Relation  (1.22)  and  the  embeddings  MWp(Q,d)  t— >■  MWp(0,d)  t— >■  MWp(0,d) 
yield  the  statement  of  Theorem  1.1. 


10 


Remark  1.2.  Let  1  <  p  <  oo  and  r  >  1/p.  Then 

sm(Mw; )  x 

The  upper  estimate  follows  from  (1.16).  The  lower  estimate  follows  from  Theo¬ 
rem  1.1. 

2.  Optimal  cubature  formulas  and  nonlinear  approximation 

The  relations  (1.10)  and  (1.11)  can  be  interpreted  as  a  connection  between  the 
error  of  the  cubature  formula  (A,  £)  on  the  class  MW J  and  the  approximation  error 
of  a  special  function  1  =  fTd  Fr(x)d/j,  by  to- term  linear  combination  of  functions 
Fr(£i  —  •),  j  =  1, . . . ,  to.  The  latter  problem  is  a  problem  of  nonlinear  to- term 
approximation  with  regard  to  a  given  system  of  functions,  in  the  above  case  with 
regard  to  the  system  {Fr(x  —  -),x  G  Td}.  The  problem  of  nonlinear  m-term  ap¬ 
proximation  has  attracted  a  lot  of  attention  during  last  ten  years  because  of  its 
importance  in  numerical  applications  (see  surveys  [15]  and  [51]).  In  this  section 
we  will  use  some  known  results  from  to- term  approximation  in  Banach  spaces  for 
estimating  the  error  of  optimal  cubature  formulas.  Let  1  <  q  <  oo.  We  define  a 
set  JCq  of  kernels  possessing  the  following  properties.  Let  and  f ly  be  measurable 
sets  for  variables  x  and  y  respectively.  Let  K(x,y )  be  a  measurable  function  on 
x  Lly.  We  assume  that  for  any  i  6  flj.  K{x ,  •)  G  Lq(Lly),  for  any  y  G  Vty  the 
K(-,y )  is  integrable  over  and  fn  K(x ,  • )dx  G  Lq(Lly).  For  a  kernel  K  G  /Cp/  we 
define  the  class 

Wp  :=  {/  :  /  =  K(x,y)cp{y)dy,  M\Lp(ny)  <  !}•  (2.1) 

Then  each  /  G  is  integrable  on  Llx  (by  Fubini’s  theorem)  and  defined  at  each 
point  of  Ltx.  We  denote  for  convenience 

J{y)  ■■=  Jk(v)  ■=  /  K(x,  y)dx. 

For  a  cubature  formula  Am(-,£)  we  have 

r,  m 

A  m{Wp,t)=  sup  \  (J{y)  -  ^»K(^,y))(p{y)dy\  = 

H^lUp(Oy)<l  d  £}y  ^=1 

m 

=  l|.?(0  -  E  vne. Olivia,).  (2.2) 

n=l 

We  use  the  same  as  above  definition  of  the  error  of  optimal  cubature  formula  with 
to  knots  for  a  class  W 


Sm(W)  : 


inf  A  m(W,0- 


Thus  by  (2.2) 


U<)=  inf  IU(-)-EA(‘X«'‘'-)llv«l»)- 

..  —  1 

e,..,r  * 

We  will  now  introduce  some  notations  and  concepts  from  the  theory  of  to- term 
approximation  in  Banach  spaces. 

Let  Xbea  Banach  space  with  norm  1 1  •  1 1 .  We  say  that  a  set  of  elements  (functions) 
V  from  X  is  a  dictionary  if  each  g  G  D  has  norm  one  (||g||  =  1), 

g  EV  implies  —  g  G  T>, 


and  span'D  =  X. 

We  will  discuss  in  this  section  two  types  of  greedy  algorithms  with  regard  to  D. 
For  an  element  /  6  X  we  denote  by  Nf  a  norming  (peack)  functional  for  /: 

11^/11  =  1.  Nf(f)  =  11/11- 

The  existence  of  such  a  functional  is  guaranteed  by  Hahn-Banach  theorem.  Let 

t  '■=  {tfcjfcli  be  a  given  sequence  of  positive  numbers  tk  <  1,  k  =  1, _ We  define 

first  the  Weak  Chebyshev  Greedy  Algorithm  (WCGA)  that  is  a  generalization  for 
Banach  spaces  of  Weak  Orthogonal  Greedy  Algorithm  defined  and  studied  in  [49] 
(see  also  [16]  for  Orthogonal  Greedy  Algorithm). 

Weak  Chebyshev  Greedy  Algorithm  (WCGA).  We  define  :=  f^T  :=  f. 

Then  for  each  to  >  1  we  inductively  define 
!)■  Vcm  :=  V™  e  V  is  any  satisfying 

Nfcm_t  >  tm  SUP  Nfc .  (g). 

g£V 

2) .  Define 

:=  :=  span{^}™=1, 

and  define  G cm  :=  G\ to  be  the  best  approximant  to  /  from  $m. 

3) .  Denote 

i*C  . _  fC,T  . _  T  _  C(C 

J  m  *  Jm  *  J  ^ m ' 

We  define  now  the  generalization  for  Banach  spaces  of  the  Weak  Relaxed  Greedy 
Algorithm  studied  in  [49]  in  the  case  of  Hilbert  space. 

Weak  Relaxed  Greedy  Algorithm  (WRGA).  We  define  /q  :=  /q’t  :=  /  and 
G o  :=  GqT  :=  0.  Then  for  each  m  >  1  we  inductively  define 
!)■  :=  eV  is  any  satisfying 

N (Pin  ~  l)  >  tm  SUp  Nf^  (g  -  G^_j). 
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2) .  Find  0  <  Am  <  1  such  that 

11/  -  ((1  -  \m)Grm-i  +  Am^)ll  =  inf  11/  -  ((1  -  a )Grm_1  +  A^)|| 

U  \  A  \  1 

and  define 

Grm  :=  G%  :=  (1  -  A m)Grm_1  +  \m<prm. 

3) .  Denote 

/r  . _  /?r*,T  . _  n  s~ir 

m  *  J  m  *  J  ^  m ' 

The  term  ’’weak”  in  both  definitions  means  that  at  the  step  1).  we  do  not  shoot  for 
the  optimal  element  of  the  dictionary  which  realizes  the  corresponding  sup  but  are 
satisfied  with  weaker  property  than  being  optimal.  The  obvious  reason  for  this  is 
that  we  don’t  know  in  general  that  the  optimal  one  exists.  Another  practical  reason 
is  that  the  weaker  the  assumption  the  easier  to  satisfy  it  and  therefore  easier  to 
realize  in  practice. 

We  present  in  this  section  results  of  convergence  and  the  rate  of  convergence  for 
the  two  above  defined  methods  of  approximation.  It  is  clear  that  in  the  case  of 
WRGA  the  assumption  that  /  belongs  to  the  closure  of  convex  hull  of  V  is  natural. 
We  denote  the  closure  of  convex  hull  of  V  by  AifiD).  It  has  been  proven  in  [49] 
that  in  the  case  of  Hilbert  space  both  algorithms  WCGA  and  WRGA  give  the 
approximation  error  for  the  class  A\(T>)  of  the  order 

m 

(i  +  X^r1/2- 

k= i 

We  discuss  here  approximation  in  uniformly  smooth  Banach  spaces.  For  a  Banach 
space  X  we  define  the  modulus  of  smoothness 

p{u)  :=  sup  {h\\x  +  uy\\  +  \\x-uy\\)-l). 

IMMMN1  2 

The  uniformly  smooth  Banach  space  is  the  one  with  the  property 


lim  p(u)/u  =  0. 

it— >-0 

The  following  convergence  result  has  been  proven  in  [50] .  In  the  formulation  of  this 
result  we  need  a  special  sequence  which  is  defined  for  a  given  modulus  of  smoothness 
p(u)  and  a  given  r  =  {tk}^=1- 

Definition  2.1.  Let  p(u )  be  an  even  convex  function  on  (—00,00)  with  the  prop¬ 
erty:  p( 2)  >  1  and 

lim  p(u)/u  =  0. 

u—¥  0 

For  any  r  =  {tk}<£=1,  0  <  t^  <  1,  and  0  <6  <  1/2  we  define  :=  £ rn(p,r,6 )  as  a 

number  u  satisfying  the  equation 


p(u)  =  Otm.U. 


(2.3) 
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Theorem  2.1.  Let  X  be  a  uniformly  smooth  Banach  space  with  the  modulus  of 
smoothness  p(u).  Assume  that  a  sequence  r  :=  satisfies  the  condition:  for 

any  6  >  0  we  have 

OO 

E  tm£,m{p,T,6)  =  OO. 
m=  1 

Then  for  any  f  e  X  we  have 


lim  ll/mTH=0. 

ra— ^  oo 

Corollary  2.1.  Let  a  Banach  space  X  have  modulus  of  smoothness  p(u)  of  power 
type  1  <  q  <  2;  (p(«)  <  qu9).  Assume  that 

OO 

E  C  =  °°,  P  =  (2-4) 

m=l  y 

Then  WCGA  converges  for  any  f  &  X. 

It  is  well  known  (see  for  instance  [14],  Lemma  B.l)  that  in  the  case  X  =  Lp, 
1  <  p  <  oo  we  have 


(u)  <  i  UP^P  if  1  -  p  -  2’ 

P  ~  \  (p  —  l)u2 /2  if  2  <  p  <  oo. 

It  is  also  known  (see  [28],  p.63)  that  for  any  X  with  dim  X  =  oo  one  has 

p(u)  >  (1  +  u2)1/2  -  1 


and  for  every  X ,  dim  X  >  2, 


p(u)  >  Cu 2,  C  >  0. 

This  limits  power  type  modulus  of  smoothness  of  nontrivial  Banach  spaces  to  the 
case  1  <  q  <  2. 

Let  us  apply  Corollary  2.1  for  numerical  integration.  Consider  a  dictionary 

V  :=  V(K, p')  :=  {p  :  g(x,y )  =  ^(®,J/)/||A'(a;,-)IUJ),(n1/)} 

(in  case  ||/C(x,  -)I|l  /(n„)  =  0  we  set  g(x,  •)  =  0),  and  define  a  Banach  space  X  := 
X(K,p')  as  the  Lp>  (Qy) -closure  of  span  of  D.  Assume  now  that  Jk  €  X.  Then  for 
1  <  p'  <  oo  the  WCGA  satisfying  (2.4)  with  q  =  min(2,p/)  provides  a  deterministic 
algorithm  of  constructing  a  sequence  of  cubature  formulas  A^(-,£)  such  that 

^rn^Wp  , 0  0  as  m—>oo. 

We  will  discuss  in  more  detail  a  question  of  the  rate  of  convergence.  The  following 
theorem  has  been  proven  in  [50]. 
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Theorem  2.2.  Let  X  be  a  uniformly  smooth  Banach  space  with  the  modulus  of 
smoothness  p(u )  <  ~fuq ,  1  <  q  <  2.  Then  for  a  sequence  r  :=  {tk}kLi,  tk  G  [0,1], 
k  =  1,2,...,  we  have  for  any  f  G  ^4i(T>)  that 

m 

\uYi\<c,tq,~,)(i  +  Y.tir1/r, 

k=  1 


m 

ll/ni<C'2(9,7)(l+E‘S)“1/’’-  p  :=  “ TT' 

fc=l  ^ 

with  constants  Ci(q,  7),  «  =  1,  2,  which  may  depend  only  on  q  and  7. 

Corollary  2.2.  In  a  particular  case  r  =  tk  =  t,  k  =  1,2, ,  with  some 

t  G  (0, 1],  we  have  under  assumptions  of  Theorem  2.2  that 

ll/£TH<Ci(<;,7,<)m-1/p, 

ll/rnT|l  <  p  := 

In  order  to  apply  Theorem  2.2  in  numerical  integration  for  IL^f  we  need  to 
check  that  Jk  G  *4.i(X,(-£f,p/))  (or  there  exists  a  positive  constant  c  such  that 
cJx  G  Ai(T>(K,p')).  It  could  be  a  difficult  problem.  An  inspection  of  the  proof  of 
Theorem  2.2  shows  that  it  is  sufficient  to  check  that 


dx  <  00. 


We  formulate  this  as  a  theorem. 

Theorem  2.3.  Let  W ^  be  a  class  of  functions  defined  by  (2.1).  Assume  that 
K  G  Kip'  satisfies  the  condition 


\\K(x,-)\\Lp,{ny)dx  <  M 


and  Jk  G  X(K,p').  Then  for  any  m  there  exists  (provided  by  WRGA  with  r  =  {t} ) 
a  cubature  formula  Am(-,£)  with 


E  ivi  <  m 

u=l 


and 


a  m(w«o<MC(P,nx,ny,t) 


m 


m 


1/2 

) 

-1  /p 


1<P<2, 

2  <  p  <  00. 


Let  us  consider  a  particular  example  of  K(x,  y )  =  (2n)  dF(x—y),  Llx  =  Qy  =  Td. 


We  denote  the  corresponding  class  Wff  by  W 
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Proposition  2.1.  Let  1  <  p  <  oo  and  ||F||p/  <  M.  Then  the  kernel  K(x,y)  = 
(2tt )~dF(x  —  y)  satisfies  the  assumptions  of  Theorem  2.3. 

Proof.  It  is  obvious  that  K  e  JCP>.  Next, 

{2ir)-d  [  || F(x  -  -)|| P'dx  =  \\F\\p>  <  M. 

Jjd 

It  remains  to  check  that  Jk  £  X(K,p').  We  have 


Jk(v)  =  (27 r) 


ljd 


F(x  —  y)dx  =  F(  0). 


Denote 


Sn{F,x):=  F(k)ei<yk,x\ 

I  kj\<N,j=l,...,d 

Then  it  is  well  known  (by  the  M.  Riesz  theorem)  that 

\\F  —  Sat (P)  | |p'  — >■  0  as  N  — >  oo. 

For  a  given  N  we  consider  a  cubature  formula 

N  N 

qN(f)  ■■=  N-d  E-E  f(2'K/jll/N,...,2'Kiad/N). 

Ud  =  1  Mi  =  l 

Then  we  have 

\\F(0)-qN{F(.-y))\\p,= 

=  ||-F(0)  -  qN(SN^(F,  -~y))+  qN(SN^(F,  .-y)-Ffi-  y))  ||p,  = 

=  \\qn(Sn-i(F,  ■  -y)  -  F(-  -  y)) ||p/  <  ||Sjv-i(-F)  -  F||p/  ->■  0 

as  N  — >  oo.  This  proves  the  proposition. 

Theorem  2.3  and  Proposition  2.1  yield  the  following  result. 

Theorem  2.4.  Let  1  <  p  <  oo.  Assume  F  e  Lp>(Td),  p1  =  p/{p  —  1).  Consider 
the  class 

WF:={f:f  =  F*cp,  |M|P<1}. 


Then  for  any  m  there  exists  a  cubature  formula  Am(-,£)  with  X^!T=i  I'Vl  —  11-^1 


and 


Am(W'pF,0<C(p,d)||F||, 


m 


m 


1/2 

) 

-1/p 


1<P<2, 

2  <  p  <  oo. 


A  sequence  of  {Am(-,  £)}  from  Theorem  2.4  can  be  obtained  by  applying  WRGA 
with  a  fixed  r  =  {/} .  We  will  describe  this  procedure  in  detail.  Denote  q  :=  p1 .  We 
have  1  <  q  <  oo.  Then 


X (K,  q)  =  span{T(a;  —  0,  %  £  Td}, 
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with  the  closure  in  Lq(Td).  It  is  well  known  that 

span {F(x  -  -),x  G  Td}  =  L^T1)  n  spanje^),  k  :  F(k )  /  0}. 

As  a  dictionary  we  have 

V  =  {±F(x--)/\\F\\p,,xeTd}. 

We  note  that  in  the  proof  of  Proposition  2.1  we  actually  proved  that 

FW/HFIU  6^(15). 

Let  us  for  simplicity  give  an  algorithm  for  F  satisfying  ||_F||p/  =  1  (otherwise  we 
take  F/\\F\\pt).  We  begin  with  /  =  F( 0).  First,  we  construct  a  norming  functional 
Nf.  It  is  known  that  for  1  <  q  <  oo  the  Nf  acts  as 

Nf(g)  =  (2n)-d  [  ||/|| J  ^ |/|9  2/^.  (2.5) 

JTd 

Fix  t  G  (0, 1)  and  find  satisfying 

\Nf{F(£1  —  -))|  >  t  sup  \  Nf(F(x  —  •))!• 

X 

Find  Ai  G  [0, 1]  and  =  ±1  such  that 

||F(0)  -  \1e1F(t1  -•)!!,=  |nf  HF'(O)  -  A F^1  -  •)||(?. 

Denote 

G1:=X1e1F(e~-)  and  /,  :=  F(0)  -  Gx. 

We  now  describe  the  mth  step.  Assume  we  have  already  obtained  Grn  \  and  frn  l  • 
Define  by  (2.5)  the  norming  functional  Nfm_1  anf  find  satisfying 

I Nu-AFir  -  -))l  >  *sup  | Nf^iFix  -  0)|. 

X 

Denote  em  =  sign Nfm_1(F(£m  —  •))  and  find  A m  G  [0, 1]  such  that 

\\F(0)  -  ((1  -  A m)Gm_!  +  A memF(r  -  -))||g  = 

=  inf  ||-F(0)  —  ((1  —  \)Gm-i  +  AemF'(^m  —  -))\\q. 

Ag[0,1] 

Denote 


Gm  (1  —  \m)Gm-i  +  AmemF1(^m  —  -)j  fm  f  —  Gr 


After  to  steps  we  obtain  a  cubature  formula  with  knots  £  =  (£-*■, . . . ,  £m)  and  weights 
A  =  (eiAi, . . . ,  emAm).  Nice  properties  of  this  construction  are:  (1)  -V  —  lj 

(2)  proceeding  from  the  step  to  —  1  to  the  step  to  we  add  one  new  knot  and 
change  in  a  simple  way  weights  e^A k  from  the  previous  step. 
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3.  Lower  estimates  for  the  classes  MW£ 


We  will  present  here  some  methods  of  obtaining  lower  estimates  of  hm(W,  £)  for 
the  classes  MW 

It  will  be  convenient  for  us  to  assume  that  the  functions  are  1-periodic  and  to 
keep  the  notation  f Id  =  [0,  l]d.  We  begin  with  the  following  result. 

Theorem  3.1.  For  any  cubature  formula  (A,  £)  with  m  knots  the  following  relation 
holds  (r  >  1/2) 

A m(MWT,0  >  C'(r,d)m-r(logm)V. 


Proof.  Let  us  denote 

m 

A  (A:)  =  Am(ei27r(fc’a:),0  =  J2 

n=i 

for  the  cubature  formula  (A,  £)  and  k  E  Zd. 

Lemma  3.1.  The  following  inequality  is  valid  for  any  r  >  1/2 

A(k) \2i,(k)-2r  >  C(r,  d)  |  A(0)  |2m_2r’(logm)d_1, 

k^O 

where  kj  :=  max(|/cj|,  1)  and  u(k)  =  ]Xj=i  kj  ■ 

First  we  will  deduce  Theorem  3.1  from  Lemma  3.1  and  then  prove  this  lemma. 
We  assume  that  |A(0)|  >  1/2  because  otherwise  it  is  sufficient  to  take  as  an 
example  the  function  f(x)  =  1.  Let  us  consider  the  function 

f(x)  =  ^  A {k)iy(k)-2rei2lv^x) 

k^O 


where  A (k)  is  the  complex  conjugate  to  the  A (k).  Then 

ii/wii2= 

\W0  J 

and 

a  m(/,o  -  /( o)  =  £|a(*o|  Wr2r- 

k^  0 

By  (3.1)  and  (3.2) 

a m(Mwz,o>  (5>(£o|2i#r2r  I 
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(3.1) 

(3.2) 


using  Lemma  3.1  we  get 


A m(MW;,0  >  C(r,d)m-r(logm)^, 
which  proves  Theorem  3.1. 

We  turn  to  the  proof  of  Lemma  3.1.  Let  b(x )  be  an  infinitly  differentiable  function 
such  that  b(x)  =  0  for  x  ^  (0, 1)  and  b(x)  >  0  for  x  e  (0, 1).  Let  m  be  given,  choose 
n  e  N  such  that 

2m  <  T  <  4m  (3.3) 

Denote  for  s  =  (si, . . . ,  sd ),  Sj  nonnegative  integers, 

d 

b,(x)  =  Y[b(2‘‘+2xj), 

3  =  1 


and 

Ys  =  {y  e  Qd  such  that  A m(bs(x  -y),£)  =  0}. 

It  is  easy  to  verify  that  for  all  s  with  ||s||i  =  n,  the  estimate 

|W|>C'(d)>0 

is  valid  for  the  measure  |TS|  of  the  set  Ys.  Further, 

|6S(0)A(0)|2|W|  =  [  \Am(bs(x-y),Z) -bs(0)A(0)\2dy< 

Jy3 

<  f  \Am(bs(x  -  y),£)  -bs(0)A(0)\2dy  =  ^2\A(k)\~\bs(k)\2 .  (3.4) 

J^d  k^0 

Let  a  :=  [r]  +  1.  Then  for  s  such  that  ||s||i  =  n  we  have 

d 

|M*0|  <  C(d,  a)  ]^[  (2~Sj  min(l,  2aSj  (kj)~a)  = 

3=  1 

d 

=  C{d,  a) 2(r-Vn  JJ  2~rSj  min(l,  2as>  (kj)~a) .  (3.5) 

3= 1 

By  summing  the  relation  (3.5)  over  all  s  such  that  ||s||i  =  n  and  using  the  inequal¬ 
ities 

d 

II  2~rSj  TSj  fa)  “)  < 

||4*=nJ  =  l 

d  oo  d 

<I1E  2~"‘  ndn(l,2"ifo)-“)  «  n<b)- 

3  1  «./  -0  j  1 
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we  get  from  (3.4) 


nd-12-2"|A(0)|2C'(d)  <  22(r-VnC(d,r)  ^|A(fc)|V(/c)-2r, 

k^O 


which  proves  the  lemma. 

We  proceed  to  a  stronger  result  than  Theorem  3.1.  We  will  replace  the  class 
MWf  in  Theorem  3.1  by  MW, J  with  2  <  p  <  00.  We  prove  the  following  proposi¬ 
tion. 

Theorem  3.2.  The  following  lower  estimate  is  valid  for  any  cubature  formula 
(A,  £)  with  to  knots  (r  >  1/p) 

Am(MWp,£)  >  C(r,  d,p)m~r (\ogm)  — ,  1  <  p  <  00. 


Proof.  The  proof  of  this  theorem  is  based  on  the  following  theorem  on  existence  of 
generalized  Rademacher  type  polynomials.  We  introduce  first  some  notation.  Let 

n(N,d)  =  {(ai, . . .  ,ad)  G  Rd  :  <  Nj,  j  =  l,...,d}, 

where  Nj  are  nonnegative  integers  and  N  :=  (N\, . . . ,  Nd).  We  denote 

T(N,  d)  :=  {t  :t=  J]  ckei{k'x)}. 

fcen(N,<i) 


Then 

d 

dimT(N,  d)  =  JJ(2A^-  +  1)  =:  tf(N). 

3= 1 

Theorem  3.3.  Let  e  >  0  and  a  subspace  T  C  T(N,  d)  be  such  that  dim  T  >  e?9(N). 
Then  there  is  a  t  G  T  such  that 


Halloo  =  1,  ll^lb  >  C(e,d)  >  0. 

The  proof  of  this  theorem  is  based  on  the  lower  estimates  of  the  volumes  of  the 
sets  of  Fourier  coefficients  of  bounded  trigonometric  polynomials  from  T(N,  d)  (see 
[45]  and  [47]). 

First,  we  prove  the  following  assertion. 

Lemma  3.2.  Let  the  coordinates  of  the  vector  s  be  natural  numbers  and  ||s||i  =  n. 
Then  for  any  N  <  2n~1  and  an  arbitrary  cubature  formula  (A,  £)  with  N  knots 
there  is  a  ts  G  T((  2Sl , ...  ,2  Sd),d)  such  that  ||ts||oo  <  1  and 

ts(0)-AN(ts,O>C(d)>0. 
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(3.6) 


Proof  of  Lemma  3.2.  Let  N  and  (A,  £)  be  given.  Let  us  consider  in  T((2Sl  1,...,2Sd  1),d) 
the  linear  subspace  T  of  polynomials  t  satisfying  the  conditions 

t(aj)=  o,  j  =  l  (3.7) 


Then 


dim  T  >  2 n  -N>  2n~1. 


Consequently,  by  Theorem  3.3  there  is  t1  G  T  such  that  ||i1||00 

H^lb  >  C(d)  >0. 


1  and 


(3.8) 


We  difine 

^  Ef=i^<l/2, 

ts(x)  —  <  o 

(  |t1(s)|  otherwise. 

The  relations  (3.7)  and  (3.8)  prove  the  lemma. 

We  now  complete  the  proof  of  Theorem  3.2.  Let  to  be  given.  We  choose  n  such 
that 

to  <  2n~1  <  2 TO. 

We  consider  the  polynomial 

t(x)  =  E  ts(x^ 

llslh=™ 

where  ts  are  polynomials  from  Lemma  3.2  with  N  =  m.  Then 

i(0)-Am(t,0>C(d)nd-1.  (3.9) 

Let  us  estimate  ||t^||p,  2  <  p  <  oo.  We  will  use  the  following  corollary  of  the 
Littlewood-Paley  theorem.  Let  f  E  L i;  denote  for  jieZj 

p(u)  :={k:k  =  (k1,...,kd),  [2Uj~1]  <  \k3 1  <  2u\j  =  1, . . . ,  d},  (3.10) 

Su(f):=  E  f^V{k’X)- 

k£p(u) 

Then  for  /  E  Lp,  2  <  p  <  oo  one  has  the  inequality 

ll/ll,  <  CM (£  |K(/)I|J)1/2. 


fW||  < 

P 


£  IM«<r)) 

llull  i<ri-\-d 
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We  have 


1/2 


Using  the  Bernstein  inequality  we  continue  the  estimate 

1/2 

<1  22r^u^1  ||5„(£)| 

\IUIIi<n+(i 

Next  we  have 

Su(t)  =  ^2  MU)  =  E  mo- 


2 

P 


s  i=n 


llsl|i=™ 

S  +  1>M 


By  the  inequality  ||ts||oo  <  1  we  get  from  here 

||Mt)||p  <  (n  +  d  +  1-  ||u||i) 
The  estimates  (3.11),  (3.12)  result  in 


d- 1 


1/2 


|*(r)||p«  (  E  22HMU(n  +  d+l-H1) 

||it||i  <n+d 


2(d-l) 


(3.11) 


(3.12) 


— .  (3.13) 


Comparing  (3.9)  and  (3.13)  we  get  the  conclusion  of  Theorem  3.3  for  2  <  p  <  oo. 
Clearly,  the  lower  estimate  for  1  <  p  <  2  follows  from  the  estimate  which  we  have 
just  proved. 

Theorem  3.2  gives  the  same  lower  estimate  for  different  parameters  1  <  p  <  oo. 
It  is  clear  that  the  bigger  the  p  the  stronger  the  statement.  We  now  discuss  an 
improvement  of  Theorem  3.2  in  the  particular  casep  =  1.  We  will  improve  the  lower 
estimate  by  replacing  the  exponent  (d—  l)/2  by  d—  1.  However,  this  improvement 
will  be  proved  under  some  (mild)  assumptions  on  the  weights  of  a  cubature  formula 
(A,  £)  and  also  for  a  slight  modification  of  the  classes  MW{.  Denote 


Frfi(x)  :=  1  +  2  ^ 2  k  r  cos  kx,  iel; 

k= 1 


d 

Fr, o(x)  ■=  JJ  Frfi(xj),  x  =  (a;i, . . . ,  xd)  G  Td; 
j= i 

MWl„  :={/:/  =  f  *  Fr,0,  |M|i  <  1}. 

It  is  clear  that  in  the  case  r  is  an  even  integer  we  have  MW[Q  =  MW[.  Let  B  be 
a  positive  number  and  Q(B,m)  be  the  set  of  cubature  formulas  Am(-,£)  satisfying 
the  additional  condition 

m 

ElA«lss. 

/'  A 

We  will  obtain  the  lower  estimates  for  the  quantities 


C(W)  := 


inf 


Am(W,0- 


We  will  prove  the  following  relation. 
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Theorem  3.4.  Let  r  >  1.  Then 


C.{MWl0)  >  C{r,B,d)m-r(\ogm)d-\  C(r,B,d)  >  0. 

Proof.  We  use  a  similar  to  the  above  notation 

m 

A(k)=Am(ei^x\0  =  J2X^{k,n- 

n=i 

In  the  case  |A(0)|  <  1/2  it  is  sufficient  to  consider  a  function  f(x)  =  1  as  an 
example,  and  therefore  we  will  assume  that  |A(0)|  >  1/2.  Considering  the  cubature 
formula  Am(-,£)  =  Am(-,  £)/A(0),  we  see  that 

Am(W,0  >  Am(W,0/4 

for  W  such  that  ^(/  —  /( 0))  E  W  provided  f  E  W  and  A m(-,£)  is  exact  on  the 
function  f(x)  =  1,  i.e.  A(0)  =  1.  Thus  it  is  sufficient  for  our  purpose  to  consider 
the  cubature  formulas  Am(-,£)  satisfying  the  additional  condition  A(0)  =  1.  Let  us 
consider  the  cubature  formula  A'  constructed  with  the  use  of  Am(-,£)  as  follows: 

m 

L>~\ 

where  Xu  is  the  complex  conjugate  to  the  \v.  Then 

m  m 

A'(Jfe)  :=  A' (ei('k,x'>)  =  =  \A{k)\2 .  (3.15) 

v=\  At=l 

The  function  Fr  {)  belongs  to  the  closure  in  the  uniform  norm  (r  >  1)  of  the  class 
MW{ o-  Consequently,  by  (3.15)  and  Lemma  3.1,  we  obtain 

A'(Mlb/0)  >  A'(Fr,0)  -  Fr, 0(0)  =  ^  A'(k)Fr,0(k)  = 

k^O 

\  A{k)\2 Fr^0(k)  >  C{r,d)m-r (log m)d_1.  (3.16) 

k^O 

On  the  other  hand,  for  the  cubature  formula  A'  we  have 

m 

A '(/)  -  /(0)  =  V  A„(A m(f(x  -  0,0  -  /(0)) 

V=1 

which,  for  Am(-,£)  E  Q(B,N),  implies  the  inequality 

A/iffif^)  <  BAm(MWl O,0-  (3-17) 

Relations  (3.16)  and  (3.17)  yield  the  required  lower  estimate  for  (5^(MWf  0). 

Let  us  discuss  how  Theorems  3.2  and  3.4  can  be  used  for  estimating  from  below 
the  generalized  discrepancy  Dr(£,  A,  m,  d)q.  Theorem  3.2  combined  with  Theorem 
1.1  and  Remark  1.2  implies  the  following  result. 
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Theorem  3.5.  Let  1  <  q  <  oo  and  r  be  a  positive  integer.  Then  for  any  points 
£  =  (£1, . . .  ,  £m)  C  Lid  and  any  weights  A  =  (Ai, . . . ,  Am)  we  have 

Dr(£,  A,  to,  d)q  >  C(d,  r)mAr  (logm)^1^2 

with  a  positive  constant  C(d,r). 

We  now  turn  to  application  of  Theorem  3.4.  Let  r  be  an  even  integer,  then 
MWf0  =  MWf.  Assume  that  the  given  cubature  formula  Am(-,£)  G  Q(B,m). 
Then  using  the  definition  of  _Dr(£,  A,  to,  d)oo  (see  (1.15))  and  the  embedding 

MWp(Lld)  t— >■  MWf(nd)  we  get 

Dr(£,A,m,d)oo  »  Am(MW'i(fid))0-  (3-18) 

By  (1.22)  and  the  embedding  MWf  c— >■  MWf  (XI  d)  we  obtain 

Am(M^(Od),0  »  A^(MW[,0),  (3.19) 

where  0  =  (01, . . . ,  #m),  =  — 7r  +  27177^, 

rf  =  j  =  l,...,d; 

d 

A'=A,n^;),  /i  =  1, . . .  ,  m. 
i=i 

Next,  it  is  clear  that  Am(-,£)  G  Q(B,m)  implies  that  A'm(-,6 )  G  Q(C(d)B,m). 
Therefore,  by  Theorem  3.4  we  get 

A^(MW1r,0)  >  m-^logm)^1.  (3.20) 

Combining  (3.18)  (3.20)  we  obtain  the  following  statement. 

Theorem  3.6.  Let  B  be  a  positive  number.  For  any  points  £1, . . .  ,£m  C  and 
any  weights  A  =  (Ai, . . . ,  Am)  satisfying  the  condition 

m 

E 

n=i 


we  have  for  even  integers  r 

Dr(£,  A,  m,  d)oo  >  C(d,  B,  r)m~r  (log  m)d~1 
with  a  positive  constant  C(d,B,r). 
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Corollary  3.1.  Let  r  be  an  even  integer.  Then  we  have  for  the  r -discrepancy 

Dr(£,  m,  d)oo  :=  Dr(£,  (1/m, . . . ,  1/m),  d.m)^  >  m~r(log  m)d~1. 

The  case  p  =  oo  is  excluded  in  Theorem  3.2.  There  is  no  nontrivial  general  lower 
estimates  in  this  case.  We  will  give  one  conditional  result  in  this  direction. 

Theorem  3.7.  Let  the  cubature  formula  (A,  £)  be  such  that  the  inequality 

A m(MWp,0  «:  m~r(logm)(d~1)/2,  r  >  1/p , 

holds  for  some  1  <  p  <  oo. 

Then 

Am(MW^,0  ^m-’Xlogm)^-1)/2. 

Proof.  We  denote  as  above 

m 

A(fe)  = 

l=i 

Let  us  consider  the  function 

0A,€,r(aO  =  E  A(/c)Tr(fc)e^  -  1. 
k 

Then  for  the  quantity  A we  have 
Am(MW;,0=  sup  |Am(/,0-/(0)|  = 

f£MW£ 

=  sup  |  Am  (Fr(x)  *  (p(x),£)  —  (p(0)  |  = 
llvllp^i 

=  sup  \{g\£,r(-y),v(y))\  =  ||^A,€,r||p',  P'  =P/(p-l)- 

(3.18) 

Consequently,  by  hypothesis  of  Theorem  3.7,  for  some  1  <  q  <  oo,  (g  =  p')  we 
have 

|  |t?A,£,r  1 1  q  <  m~r(logm)(<i~1)/2.  (3.19) 

Further,  for  arbitrary  1  <  a  <  b  and  /  £  fj  the  following  inequality  holds 

ll/ll.  <ii/iuii/iii-.  K  =  (y  -  y)  f1  -  y)  •  (3-20> 

By  Theorem  3.2  we  have  for  any  1  <  z  <  oo 

A m(MWJia,0  »  m~r(log m)^d_1^2.  (3.21) 

Therefore,  by  (3.18) 

||5A,€,r|U'  »  m~r(logm)(<i~1)/2.  (3.22) 

Setting  now  b  =  q,  a  =  |(6  +  1),  z*  =  a  we  get  from  relations  (3.20),  (3.19),  (3.22) 

||0A,€,r||i  >  m~r(logm)(d~1)/2. 

It  suffices  to  apply  the  relation  (3.18)  to  complete  the  proof. 

Theorem  3.7  is  proved. 
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Remark  3.2.  We  have  actually  proved  the  following  inequality.  Let  1  <  pi  <  P2  < 
oo,  then  for  any  (A,  £) 

A (r  >  1/Pl). 


4.  Upper  estimates  for  the  classes  MWf 

4.1.  The  Fibonacci  cubature  formulas.  For  periodic  functions  of  two  variables 
we  consider  the  Fibonacci  cubature  formulas 

bn 

®n(f)  =  bf1  Y  f(2^^/bn,  27r{/i6n_i/6n}), 
fX=l 


where  bo  =  b\  =  1,  b„  =  bn  i  +  &n_2  are  the  Fibonacci  numbers  and  {x}  is  the 
fractional  part  of  the  number  x. 

For  a  function  class  W  we  denote 


$n(W) 


sup  I  $„(/) 
few 


(2vr)  2  J f(x)dx\, 

J2 


where  T2  =  [0,  27t] 2  is  the  period  square. 

The  following  known  result  gives  the  order  of  <hn(Wf)  for  all  parameters  1  < 
p  <  oo,  r  >  1/p. 

Theorem  4.1.  We  have 


(  bn  r (log  bny/2, 


bn  r  iog  bn, 

bnr(\og  bny-r, 

bnr(  (log  bn )  (log  log  bn ) ) x!2 , 


1  <  p  <  oo,r  >  max 

p=l,r>  1; 
2<p<oo,  i<r<|; 

2  <  p  <  oo,  r  =  1/2. 


The  lower  estimates  provided  by  Theorem  3.2  and  the  upper  estimates  from 
Theorem  4.1  show  that  the  Fibonacci  cubature  formulas  are  optimal  (in  the  sense 
of  order)  among  all  cubature  formulas  in  the  case  1  <  p  <  oo,  r  >  max(l /p,  1/2): 


6bJMW;)  x  4>„(MWg  x  6;r(l0gfc„)1''2. 


Theorem  3.4  combined  with  Theorem  4.1  implies  that  the  Fibonacci  cubature  for¬ 
mulas  are  optimal  (in  the  sense  of  order)  among  formulas  satisfying  an  additional 
restriction  <  B  in  the  case  r  an  even  integer  and  p  =  1 

S?JMW[)  x  <6„(MW[)  x  b-Tlogbn. 

For  all  other  values  of  parameters  p  and  r,  r  >  1/p ,  the  right  order  of  Sm(MWp )  is 
unknown. 
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4.2.  The  Korobov  cubature  formulas.  It  will  be  convenient  for  us  to  denote 
vectors  in  7hd  and  M.d  by  bold  letters.  Let  m  E  N,  a  =  (ai, . . .  ,ad)  E  7Ld .  We 
consider  the  cubature  formulas 


m 


fx=l 


which  will  be  called  the  Korobov  cubature  formulas. 
In  the  case  d  =  2,  m  =  bn,  a  =  (1,  b„  \ )  we  have 


Pm(/,a)  =  $„(/). 


We  note  that  in  the  case  d  >  2  the  problem  of  finding  concrete  cubature  formulas 
of  the  type  Pm(/,  a)  as  good  as  the  Fibonacci  cubature  formulas  in  the  case  d  =  2 
is  unsolved.  The  results  of  this  subsection  deal  with  the  case  d  >  2  and  are  not  as 
complete  as  the  results  of  Subsection  4.1. 

We  first  prove  an  auxiliary  assertion.  Denote 


d 

r (N)  :=  {k  =  (ku  •  •  • ,  fed)  E  :  J]  maxflfyl,  1)  <  IV}. 

i=i 

For  a  finite  set  E  the  cardinality  of  E  will  be  denoted  by  \E\. 

Lemma  4.1.  Let  n,  k,  L  be  a  prime,  a  positive  real  and  a  natural  number,  respec¬ 
tively,  such  that 

\T(L)\  <(n-l){l-2-K)/d.  (4.1) 

Then  there  is  a  natural  number  a  6  In  :=  [l,n)  such  that  for  all  m  G  T (L), 

m  0 

mi  +  am,2  4 - +  ad~1md  ^  0  (mod  n),  (4.2) 

and  the  relation  (4.2)  will  be  valid  for  all  m  E  F)(L)  :=  T(L2l)\T(L21^1),  m  /  nm', 
with  the  exception  of  no  more  than 

A\  :=  \Fl{L)\d2^K(2K  -  1)  x(ri  —  1)  x,  1  =  1,2,...  . 

Proof.  Let  a  G  /„,  be  a  natural  number.  We  consider  the  congruence 

mi  +  awi2  +  •  •  •  +  ad~1md  =  0  (mod  n).  (4.3) 

For  a  fixed  vector  m  =  (mi, . . .  ,rrid)  we  denote  by  An( m)  the  set  of  natural 
numbers  a  G  In  which  are  solutions  of  the  congruence  (4.3).  It  is  well-known  that 
for  m  ^  0,  m,j  <  n,  j  =  1, . . . ,  d  the  number  |An(m)|  of  the  elements  of  the  set 
An( m)  satisfies  the  inequality 


An  (m)  |  <  d  —  1  <  d. 
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(4.4) 


We  denote  by  G i  the  set  of  the  numbers  a  which  are  solutions  of  the  congruence 
(4.3)  for  at  least  one  of  m  G  T(L),  that  is 


G\  —  UmGr(L)An(m). 


Let  us  estimate  the  number  |Gi|  of  elements  of  the  set  G\.  By  (4.4)  and  (4.1) 
we  have 

|Gi|<  |^n(m)|  <d|r(L)|  <(n-l)(l-2-K).  (4.5) 

mer(L) 


For  any  a  G  In\G i  for  all  m  G  T(L)  we  have 

mi  +  arri2  +  •  •  •  +  ad~1md  ^  0 


(mod  n). 


Let  Gi+ 1,  l  =  1,2,...  denote  the  set  of  those  a  for  which  the  number  of  elements 
of  the  set 

Mla  :=  {m  :  m  G  F)(L),  m  /  nm7,  mi+am^H - \-ad^1md  =  0  (mod  n)} 

satisfies  the  inequality 


Then  by  (4.6) 


\K\  >  Af. 

E  Kl>^f|G1+i|. 

iGGz  +  i 


(4.6) 

(4.7) 


On  the  other  hand,  by  (4.4)  each  m  can  belong  to  at  most  d  —  1  different  sets 


MLa  and  therefore 

£  W‘J  <  d\F,(L)\. 

a£Gi  + 1 

Comparing  (4.7)  and  (4.8)  we  find 

|GI+1  <  d\F,(L)\lAi  =  (n-  1)(2“  -  l)2-'t(i+D. 


(4.8) 


(4.9) 


From  relations  (4.5)  and  (4.9)  it  follows  that 


£|G,|<n-l. 

1=1 


This  means  that  there  exists  a  number  a  G  In  which  does  not  belong  to  any  set 

Gi,  l  =  1, - This  a  is  the  required  number  by  the  definition  of  the  sets  Gi.  The 

lemma  is  proved. 

For  a  of  the  form  a  =  (1,  a, ... ,  ad _1)  denote  Pn{f,  a)  :=  Pn{f ,  a). 
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Theorem  4.2.  For  any  r  such  that  r  >  max  (1/p;  1/2)  and  any  prime  number  n 
there  is  a  natural  number  a  G  [l,n)  for  which 

Pn(MWp,a )  <  C(r,p,  d)n_r(logn)r^_1\  1  <  p  <  oo. 


Proof.  Clearly,  it  suffices  to  consider  the  case  1  <  p  <  2.  Let  k  :=  (r  —  l/p)/2,  n  be 
given  and  L  satisfies  (4.1)  and  in  addition  |T(L)|  x  n.  Let  a  be  the  number  from 
Lemma  4.1,  depending  on  n,  k,  L  and 

Zi  =  {m  :  m  G  Fi(L),  m  ^  nm',  mi+am2Ji - \-ad~1md,  =  0  (modn)}. 


Then  by  Lemma  4.1 

'  V|  Y  A  J 

and  the  error  of  the  cubature  formula  can  be  estimated  as  follows 

OO 

1=1  m£Zi 


Pn(f,a)-f( 0)  < 


I z,\  <A?, 


+ 


/(m) 


m=nm 


— :  Oy  +  <72- 


(4.10) 

(4.11) 


Let  us  estimate  ay,  02  from  (4.11)  for  /  G  MWp. 
We  denote  __ 

Mx)  =  J2  ei{k,x)- 


We  have 


o-i  = 


^1) 


1=1 


<E^«-.(/) 


pll'rdlp' ) 


Z=1 


(4.12) 


where  E^(f)p  denotes  best  approximation  of  /  in  Lp  by  trigonometric  polynomials 
with  frequencies  in  r(7V).  It  is  known  (see  [47,Ch.3,Th.3.2])  that  for  /  G  MWp, 
1  <  p  <  00  we  have 


EN{f)p  <  C(r,p,d)N~r ,  r>  0. 

Further 

n*iip'  <  <  i^r/p.  (4.13) 

From  (4.12)  by  the  estimate  of  E^(f  )p  and  by  (4.13)  and  (4.10)  we  get 

ai  r  «7i-r(logn)p(d-1>. 

For  0-2  we  have 

0-2  <  |/(nm/)|  n~r  ^^^(m/)~r|<p(nm/)|,  (4-14) 

m'^0  m' 

where  <p  is  such  that  f  =  Fr  *  ip,  ||<p||p  <  1. 

From  (4.14),  applying  the  Holder  inequality  and  the  Hausdorff-Young  theorem, 
we  get 

(\  1/p  /  \  1/p' 

^V(m/)~rpj  |^|p(nm/)|p  J  n~r||p||p  <  n~r. 

The  conclusion  of  the  theorem  follows  from  the  estimates  for  oy  and  ay. 
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4.3.  The  Frolov  cubature  formulas.  In  this  subsection  we  construct  the  opti- 

o 

mal  (in  the  sense  of  order)  cubature  formulas  for  the  classes  MWrp(Q,d ),  2  <  p  <  oo. 

The  following  lemma  plays  a  fundamental  role  in  the  construction  of  such  cuba¬ 
ture  formulas. 


Lemma  4.2.  There  exists  a  matrix  A  such  that  the  lattice  L( m)  =  Am,  where  m  is 
a  ( column )  vector  with  integer  coordinates,  has  the  following  properties 


/ 

( I'i(m)  \  \ 

L(m)  = 

: 

V 

\Ld(m) )  ) 

1° 

2° 


n-=iL;(m) 


>  1  for  all  m  7^  0; 


each  parallelepiped  P  with  volume  \P\  whose  edges  are  parallel  to  the  coor¬ 
dinate  axes  contains  no  more  than  \P\  +  1  lattice  points. 


Let  a  >  1  and  A  be  the  matrix  from  Lemma  4.2.  We  consider  the 
formula 

(A_1)Tm' 


i(a,A)(f)=  (ad\ det/l|)  1  V  / 

m£Zd 


a 


cubature 


o 

for  /  E  MW2(Lld).  Clearly,  the  number  N  of  points  of  this  cubature  formula  does 
not  exceed  C(A)ad\  det  A\. 

Theorem  4.3.  Let  a  matrix  A  be  from  Lemma  f.2  and  let  r  be  a  natural  number. 
Then 

<$>(a,A)(MWr2(Lld))  <  C  (A,  d)a-rd  (log  a)^. 


Proof.  We  will  use  the  Poisson  formula  which  we  formulate  in  a  form  convenient 
for  us.  We  denote  for  /  E  L1(Md) 


/( y)  =  [  f(x)e-2^y^dx. 

jRd 


Lemma  4.3.  Let  /(x)  be  continuous  and  have  compact  support  and  the  series 
Skezd  /(k)  converges.  Then 


5Z  -^k) =  /(n)- 

k£Zd  n£Zd 


By  Lemma  4.3  the  identity 


(ad  \  det  A|)  1 

m 


((A  ?TX)  e~27n(m’x)dx  = 


^(aAm,y)dy  =  J2f(aAm), 


(4.15) 
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m 


holds  under  the  assumption  that  the  series  in  the  right  side  of  (4.15)  converges. 
The  convergence  of  this  series  will  follow  from  further  consideration.  In  the  relati¬ 
on  (4.15)  we  carried  out  the  linear  change  of  variables  y  =  J  x. 

We  have  for  the  error  of  this  cubature  formula 

S  =  $(a,  A)(f)  -  /( 0)  =  Y  f{aAm).  (4.16) 

m^O 

We  need  the  following  simple  assertion. 

Lemma  4.4.  Let  W^pWi  <  1  and  the  support  of  ip  be  contained  in  Lid-  Then  for  any 
a  >  1  and  nonsigular  matrix  A  we  have 

y^|(p(m4m)|2  <  C(A). 


(4.17) 


(4.18) 


(4.19) 


From  the  relation  (4.17),  using  the  Cauchy  inequality  and  the  inequality  (4.18), 
we  get 

2 

|cp(a7lm)|2  =  (ad\ det  A|)  2  cm(n)  < 

n  £G 

<  (ad|  det  vl|)  2\G\  Y  |cm(n)  |"  <  (ad\ det  A\)  1C2{A)  Y,  |cm(n)|2. 

n<=G  n£G 
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Performing  the  summation  over  m  and  taking  into  account  the  relation  (4.19) 
we  get 


(f(aAm)  2  <  02(A)  (ad\  det  A|) 


(^) 


=  C2(A)  [  \ip(y)\2dy  <C2(A). 

JRd 


2 

dx  = 


The  lemma  is  proved. 

We  continue  the  proof  of  Theorem  4.3. 

0 

Let  /  e  MW 2(^0)  be  rd  times  continuously  differentiable.  We  denote  <p(x)  = 
Then  for  m  ^  0 

d 

/(aim)  =  (/(aim)  (2niaLj(m)^  r. 

3— 1 


Let  l  be  such  that 

2 i_1  <  <  2Z. 

Then  by  the  property  1°  of  Lemma  4.2  the  inequality  ||s||i  >  l  holds  for  s  such  that 
p(s)  (see  the  definition  of  p( s)  in  (3.10))  contains  a  point  aim  with  m  ^  0. 

Then 


(2vr)  rdS<  |^(aAm)|  niaL^m)l  - 

||s||i>Z  aL(m)£p(s)  j  =  l 


1/2 


^  iik-m 

\||s||i>Z  aL(m)£p(s)  j= 1 

Applying  Lemma  4.4  and  using  the  relation 

s  ni^(m)r2r«  s 2 


ip(aAm) 


1/2 


(4.20) 


2r  NT-  ^-2r||sib  (  P^ 

d 


+  1  < 


||s||i>Z  aL(m)Gp(s)  j=l  ||s||i>Z 

<  2-2rlld-1  <  a-^loga)^1. 
we  get  from  here  and  (4.20) 

S  <C  a^rd(loga)^~ . 


au 


Theorem  4.3  is  proved. 

Theorem  4.3  combined  with  Theorem  1.1  and  Theorem  3.2  implies  the  following 
theorem. 


Theorem  4.4.  Let  r  G  N  and  2  <  p  <  oo.  Then 

5m(MWp)  x  m^r  (\ogm)~^~ . 
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5.  Historical  notes,  comments,  and  some  open  problems 

First  we  will  give  a  brief  historical  survey  on  the  discrepancy.  We  refer  the  reader 
for  a  complete  survey  to  the  following  books  on  discrepancy  and  related  topics  L. 
Kuipers  and  H.  Niederreiter  [27],  J.  Beck  and  W.  Chen  [5],  J.  Matousek  [29],  and 
B.  Chazelle  [9].  We  formulate  all  results  in  the  notations  of  this  paper  and  in  the 
form  convenient  for  us.  We  use  the  following  notation 


«  d  ^  m 

D(X,  m,  d)q  '■={  |  n  - ^2x[0,a]ix^)\qda)1/q,  l<q<oo, 

j= 1  m  M=1 

d  ^  m 

D(X,  TO^oc  :=  max  117% - X  X[o,a](^M)|, 

a£  0,1  d  -LJ-  to  -<L- '  1  1 

3  =  1  fi  1 

where  X  =  (x1, . . . ,  xm ).  The  first  result  in  this  area  was  the  following  conjecture 
of  van  der  Corput  [11,12]  formulated  in  1935.  Let  E  [0, 1],  j  =  1,  2, . . . ,  then  we 
have 

lim  sup mD ((^  . . .  ,£m),rn,  ljoo  =  oo. 

ra— >-  oo 

This  conjecture  was  proved  by  van  Aardenne-Ehrenfest  [1]  in  1945: 


lim  sup 

ra— ^  oo 


log  log  log  TO 
log  log  TO 


TOT»((^1,...,^m),TO,  I)*, 


>  0. 


Let  us  denote 

D(m,  d)q  :=  inf  D(X,  to,  d)q,  1  <  q  <  oo. 

X 

In  1954  K.  Roth  [32]  proved  that 

D(m,d) 2  >  C{d)m^1{\ogm)(yd^1^2 . 


(5.1) 


In  1972  W.  Schmidt  [36]  proved 

D(m,  2)oo  >  CVnA1  logm.  (5.2) 


In  1977  W.  Schmidt  [37]  proved 

_D(to,  d)q  >  C(d,  q)m^1(log  to)^-1^2,  1  <  q  <  oo.  (5-3) 

In  1981  G.  Halasz  [21]  proved 

D(m,d) i  >  C(d)m^1  (logm)1^2 .  (5.4) 


The  following  conjecture  has  been  formulated  in  [5]  as  an  excruciatingly  difficult 
great  open  problem. 
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Conjecture  5.1  ([5]).  We  have  for  d  >  3 

D(m,d) oo  >  C(d)m^1(\ogm)d^1 . 

This  problem  is  still  open. 

We  now  present  the  results  on  the  lower  estimates  for  the  r-discrepancy.  We 
denote 

Dr(m,  d)q  :=  inf  Dr(£,  (1/m, . . . ,  1/m),  m,  d)q 

where  Dr(£,  A,m,  d)q  is  defined  in  (1.15)  and  also  denote 

D°(m,  d)q  :=  inf  Dr( £,  A,  m,  d)q. 

£  5-^- 

It  is  clear  that 

D°(m,d)q  <  Dr(m,d)q. 

The  first  result  in  estimating  the  generalized  discrepancy  was  obtained  in  1985  by 
V.A.  Bykovskii  [6] 

D°(m,d) 2  >  C(r,  d)m~r  (log  m)^1^ 2 .  (5-5) 

This  result  is  a  generalization  of  the  Roth’s  result  (5.1).  The  generalization  of  the 
Schmidt’s  result  (5.3)  was  obtained  by  the  author  in  1990,  [45],  (see  Theorem  3.5 
of  this  paper) 

D°(m,  d)q  >  C(r,  d,  q)m~r  (logm)^1^2 ,  1  <  q  <  oo.  (5-6) 

In  1994,  [48],  the  author  proved  that  for  r  even  integers  we  have  for  the  r-discrepancy 
(see  Theorem  3.6  and  Corollary  3.1  of  this  paper) 

Hr(m,d)00  >  C(r,  d)m^r  (\ogm)d^1 .  (5-7) 

This  result  encourages  us  to  formulate  the  following  generalization  of  the  Conjecture 

5.1. 

Conjecture  5.2.  For  all  d,  r  G  N  we  have 

D°(m,d) oo  >  C(r,  d)m^r  (log  m)^1 . 

The  above  lower  estimates  for  D°(m,d)q  are  formally  stronger  than  the  corre¬ 
sponding  estimates  for  D(m,d)q  because  in  D°(m,d)q  we  are  optimizing  over  the 
weights  A.  However,  the  proofs  for  D(m,  d)q  could  be  adjusted  to  give  the  estimates 
for  Di(m,  d)q.  The  results  (5.5)  (5.7)  for  the  generalized  discrepancy  were  obtained 
as  a  corollary  of  the  corresponding  results  on  cubature  formulas  (see  Theorem  1.1 
and  Theorems  3.5,  3.6).  We  do  not  know  if  existing  methods  for  D(m,  d)q  could  be 
modified  to  obtain  the  estimates  for  D°(m,d)q,  r  >  2. 

We  proceed  to  the  lower  estimates  for  the  cubature  formulas.  Theorem  3.1  and 
Lemma  3.1  were  established  in  [6].  The  proof  is  taken  from  [47].  Theorems  3.2, 
3.3,  and  Lemma  3.2  were  proved  in  [45].  Theorem  3.4  was  proved  in  [48].  Theorem 
3.7  is  from  [46].  There  are  two  big  open  problems  in  this  area.  We  formulate  them 
as  conjectures. 
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Conjecture  5.3.  For  any  d  >  2  and  any  r  >  1  we  have 

8m(MW[)  >  C(r,  d)m~r  (\ogm)d~1 . 

Conjecture  5.4.  For  any  d  >  2  and  any  r  >  0  we  have 

5m{MW^)  >  C(r,d)m-p(logm)(d-1)/2. 

We  note  that  by  Proposition  1.2,  Theorem  1.1,  and  (1.15)  Conjecture  5.3  implies 
Conjecture  5.2  and  Conjecture  5.4  implies 

D°(m,d) i  >  C(r,  d)m~r(logm)^d~1^2 .  (5-8) 

We  turn  to  the  upper  estimates.  We  begin  with  the  cubature  formulas.  We 
have  already  made  a  historical  remark  on  classes  with  bounded  mixed  derivative  in 
Section  1.  We  will  discuss  only  these  classes  here.  For  results  on  cubature  formulas 
for  the  Sobolev  type  classes  we  refer  the  reader  to  the  books  of  S.L.  Sobolev  [40],  E. 
Novak  [31],  and  the  author  [47,Ch.2].  The  first  result  in  this  direction  was  obtained 
by  N.M.  Korobov  [25]  in  1959.  He  used  the  cubature  formulas  Prn(,f,  a)  defined 
in  Subsection  4.2.  We  note  that  similar  cubature  formulas  were  also  used  by  E. 
Hlawka  [22].  The  Korobov’s  results  lead  to  the  following  estimate 

Sm(MWi)  <  C(r,  d)m~r  (\ogm)rd ,  r  >  1.  (5.9) 

In  1959  N.S.  Bakhvalov  [2]  improved  (5.9)  to 

8m{MW[)  <  C(r,  d)m~r(logm)r(d~1),  r  >  1. 

The  first  best  possible  upper  estimate  for  the  classes  MWp  was  obtained  by  N.S. 
Bakhvalov  [3]  in  1963.  He  proved  in  the  case  d  =  2  that 

6m  {MW£  )<  C(r)m-r  (log  m) 1/2 ,  re  N.  (5.10) 

N.S.  Bakhvalov  used  the  Fibonacci  cubature  formulas  defined  in  Subsection  4.1. 

In  1976  K.K.  Frolov  [18]  used  the  cubature  formulas  defined  in  Subsection  4.3 
to  extend  (5.10)  to  the  case  d  >  2  : 

(5m(MIFJ)  <  C(r,  d)m~r(logm)^d~1^2,  r  G  N.  (5-11) 

In  1985  this  estimate  was  further  generalized  by  V.A.  Bykovskii  [6]  to  r  G  M,  r  >  1. 
Bykovskii  also  used  the  Frolov  cubature  formulas.  One  can  find  these  results  in 
Section  4  of  this  paper.  Theorem  4.1  was  proved  in  [46]  and  [48].  Theorem  4.2 
and  Lemma  4.1  are  from  [44].  Theorem  4.2  addresses  the  case  of  small  smoothness: 
r  >  max(l /p,  1/2)  instead  of  r  >  1.  Lemma  4.2  is  a  well  known  result  in  algebraic 
number  theory  (see  [8]).  Theorem  4.3  was  obtained  by  K.K.  Frolov  [18,  19].  The 
proof  is  taken  from  [47].  We  note  that  there  is  no  sharp  results  for  5m(MW£)  in 
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the  case  of  small  smoothness  1/p  <  r  <  1.  It  is  an  interesting  open  problem.  The 
approach  based  on  nonlinear  to- term  approximation  (see  Section  2)  can  be  useful 
in  this  case. 

The  Frolov  cubature  formulas  [19]  give  the  following  estimate 

Sm(MWi)  <  C (r,  d)m  r  (log m)d  1 ,  r  >  1.  (5.12) 

Thus  the  lower  estimate  in  Conjecture  5.3  is  the  best  possible. 

In  1994  M.M.  Skriganov  [39]  proved  the  following  estimate 

Sm(MWp)  <  C(r,  d,p)m~r (log to) C-1)/2^  1  <  p  <  oo,  r  e  N.  (5.13) 

This  estimate  combined  with  Theorem  3.2  implies 

Sm{MW/)  ~m-r(\ogm){d-1)/2,  1  <  p  <  oo,  re  N.  (5.14) 

Another  proofs  of  (5.13)  and  Theorem  3.2  were  given  in  1995  by  V.A.  Bykovskii 

[7]- 

We  now  present  the  upper  estimates  for  the  discrepancy.  In  1956  H.  Davenport 
[13]  proved  that 

D(m,  2)2  <  Cm  ^  (log  to) . 

Another  proofs  of  this  estimate  were  later  given  by  I.V.  Vilenkin  [52],  J.H.  Halton 
and  S.K.  Zaremba  [24],  and  K.  Roth  [33].  In  1979  K.  Roth  [34]  proved 

D(m,  3)2  <  Cm-1  log  m 

and  in  1980  K.  Roth  [35]  and  K.K.  Frolov  [20]  proved 

D(m,d) 2  <  C  (d/m^1  (log  m)^1^ 2 . 

In  1980  W.  Chen  [10]  proved 

D(m,d)q  <  C(d)m^1  (logm)^-1^2 ,  q  <  00. 

The  estimate  (5.12)  and  Theorem  1.1  imply 

D°(m,d) 00  <  C(r,  d)m~r  (log  m)^1 ,  r  >  2. 

We  note  that  the  upper  estimates  for  D(m,d)q  are  stronger  than  the  same  upper 
estimates  for  D°(m,d)q. 

Let  us  also  mention  a  classical  book  of  S.M.  Nikol’skii  [30]  on  quadrature  formulas 
and  books  of  N.M.  Korobov  [26],  W.  Schmidt  [38],  and  Hua  Loo  Keng  and  Wang 
Yuan  [23]  on  discrepancy  and  related  topics.  We  discussed  in  this  paper  only  the 
case  of  the  class  MW J  of  functions  with  bounded  mixed  derivative.  There  are 
analogs  of  classes  MW//  which  are  also  natural  in  the  theory  of  cubature  formulas. 
One  can  find  results  on  numerical  integration  of  functions  with  bounded  mixed 
difference  in  the  papers  [4],  [17],  [43],  and  in  the  book  [47].  A  different  method  of 
constructing  cubature  formulas  for  functions  with  bounded  mixed  derivative  was 
suggested  by  S.A.  Smolyak  [41]  (for  further  results  see  [42]). 
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